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We study numerically the development of chimera states in networks of nonlocally coupled oscil¬ 
lators whose limit cycles emerge from a Hopf bifurcation. This dynamical system is inspired from 
population dynamics and consists of three interacting species in cyclic reactions. The complexity of 
the dynamics arises from the presence of a limit cycle and four fixed points. When the bifurcation 
parameter increases away from the Hopf bifurcation the trajectory approaches the heteroclinic in¬ 
variant manifolds of the fixed points producing spikes, followed by long resting periods. We observe 
chimera states in this spiking regime as a coexistence of coherence (synchronization) and incoher¬ 
ence (desynchronization) in a one-dimensional ring with nonlocal coupling, and demonstrate that 
their multiplicity depends both on the system and the coupling parameters. We also show that 
hierarchical (fractal) coupling topologies induce traveling multichimera states. The speed of motion 
of the coherent and incoherent parts along the ring is computed through the Fourier spectra of the 
corresponding dynamics. 
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I. INTRODUCTION 

The recent influence of the theory of networks urn on 
the classic field of coupled oscillatory units [4] has led to 
the discovery of a plethora of novel phenomena, especially 
when the coupling between units becomes more sophisti¬ 
cated, mimicking naturally interacting systems. Among 
the complex oscillatory patterns that may emerge, are 
the so-called “chimera states” EE], which have recently 
attracted a lot of attention. These are states, in which 
identically coupled units spontaneously develop coexist¬ 
ing synchronous (coherent) and asynchronous (incoher¬ 
ent) parts. 

Many recent theoretical works have focused on the 
study of chimera states in a variety of physical sys¬ 
tems. Typical models that have been numerically in¬ 
vestigated include the Kuramoto phase oscillator EE], 
periodic and chaotic maps [ini E], the Stuart-Landau 
model [13 E], the Van der Pol oscillator [12] as well as 
models addressing neuron dynamics such the FitzHugh- 
Nagumo oscillator m, the Hindmarsh-Rose model m, 
the so-called SNIPER model of excitability type-I El, 
or the Hodgkin-Huxley model [18]. Moreover, chimera 
states have been reported in populations of coupled pen- 
dula m, in autonomous Boolean networks m , in one- 
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dimensional superconducting meta-materials m, and 
time-varying networks [22] . 

Hollowing the theoretical predictions, chimera states 
were experimentally verified for the first time in popu¬ 
lations of coupled chemical oscillators [23] and in opti¬ 
cal coupled-map lattices realized by liquid-crystal light 
modulators [24]. Recently, in a purely mechanical ex¬ 
periment involving two groups of identical metronomes, 
it was shown that chimeras emerge naturally as a coex¬ 
istence of two competing synchronization patterns [25] . 
Chimeras were also realized in experiments involving 
electronic nonlinear oscillators with delay m and elec¬ 
trochemical oscillator systems [2I1128]. 

Chimera states find increasing interest due to the pos¬ 
sible connections to various phenomena observed in bio¬ 
logical and social systems. As stated in the review paper 
of Panaggio and Abrams [29] , chimera states could possi¬ 
bly explain the phenomenon of unihemispheric sleep ob¬ 
served in dolphins and some birds which sleep with half 
of their brain m as well as ventricular fibrillation m- 

In the current study, we attempt to bring in evidence 
chimera states in the fields of chemical reaction kinetics 
and ecology. Ecological models in the form of reactive 
dynamical systems have been previously used to describe 
the interactions between multiple species (referred to also 
as particles). In most cases the dynamics follows a mean 
field (ME) rate equations approach, which describe global 
dynamical features [32] [33] . In addition, numerical sim¬ 
ulations are employed to address specific details of the 
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spatial and temporal characteristics of the systems [34]- 
[36] . Along the same lines, dynamical systems describing 
interactions between different chemical species are stud¬ 
ied. Common examples involving heterogeneous catalytic 
dynamics are the reactions NO + CO [37l|38], CO + O 2 
(^[39], CO^OH [40], all taking place on a platinum cat¬ 
alyst, the reaction CO-\-2H2 on various catalysts [41], etc. 
All above processes are described as dynamical systems 
using the Langmuir-Hinshelwood mechanism [42]. Anal¬ 
ogous modeling is extensively used for the well known 
oscillatory Belousov-Zhabotinsky reaction [43| [44]. In 
epidemiology, reactive systems are frequently used where 
the species correspond to infected individuals, suscepti¬ 
ble or recovered ones, and they can dynamically infect 
other parts of the population [45H47] . All these systems 
are commonly viewed as reactive dynamics involving re¬ 
action, birth (or adsorption), and death (desorption) pro¬ 
cesses. 

In the field of experimental reactive dynamics there 
have already been several studies discussing the existence 
of chimera states. Some interesting works include the 
photosensitive Belousov-Zhabot insky reaction [23] as well 
as experiments with electrochemical oscillators [lilllH]. 
Similarly, a social analogue of a chimera state was re¬ 
ported in [ 49 ], where the authors investigated the emer¬ 
gence of localized coherence in two interacting popula¬ 
tions of social agents. Two social models were imple¬ 
mented, the Axelrod model for culture dissemination 
m and the bounded confidence model by Deffuant et 
al. [51]. Different synchronization patterns were ob¬ 
tained, including chimera states where one subpopula¬ 
tion of agents remained synchronized while the other was 
desynchronized. 

It is interesting to note here that in all experiments as 
well as in the social systems presented above, the sub¬ 
populations of elements are ab initio prepared and serve 
as a basis for the synchronized/desynchronized activity. 
On the other hand, in the numerical investigations cited 
above and in the current study, the coherent/incoherent 
regions emerge spontaneously from the interplay between 
the system parameters, dynamics, and initial conditions. 

The main ingredients for the emergence of a chimera 
state in coupled oscillatory dynamics is the stable os¬ 
cillatory pattern of the single elements together with the 
nonlocal character of the interactions. A toy model which 
exhibits well defined oscillations in the form of a limit 
cycle and is relevant in ecological, chemical, and social 
dynamics is the Lattice Limit Cycle (LLC) model, pro¬ 
posed in [48]. The LLC model is a nonlinear system that 
describes reaction, birth, and death processes and stud¬ 
ies how the concentrations of the specific species change 
with time under the influence of the control parameters. 
The LLC is an ideal model to study the development of 
chimera states in population dynamics for three reasons: 
(a) It describes common ecological/social/chemical pro¬ 
cesses using reaction, birth, and death differential equa¬ 
tions, (b) its dynamics is characterized by a limit cycle, 
and (c) when needed, it allows for direct implementa¬ 


tion on a lattice with single species occupancy for a more 
realistic representation of the dynamics. 

In this study, the capacity of the LLC model to pro¬ 
duce chimera states will be investigated. The chimera 
attributes will be discussed as the system approaches the 
critical point of the Hopf bifurcation. Furthermore, we 
will explore the influence of the connection topology on 
the chimera features. In population dynamics, exchange 
in the form of spatial diffusion between individuals resid¬ 
ing on different parts of the system is frequent: people 
move from one city to another, birds migrate, individ¬ 
ual animals circulate from one herd to another, etc. In 
natural and social dynamics, the circulation of individ¬ 
uals has often nonlocal character and thus it is natu¬ 
ral to investigate phenomena like chimera states which 
have their origin in nonlocal interactions. In particular, 
it is important to consider hierarchical connectivity in 
the coupling of the LLC dynamics, since the habitats 
of populations have often hierarchical/fractal distribu¬ 
tions and morphologies: the city locations are fractally 
distributed [52], the mammal and bird habitats are lo¬ 
cated on fragmented landscapes that affect their commu¬ 
nity dynamics [53H55] . and the viruses spread in com¬ 
plex socio-geographic networks mEIl. As we will see 
in the sequel, hierarchical connectivity induces traveling 
chimera states, a novel feature in this held. 

The work is organized as follows: In Sec. [II] we re¬ 
capitulate the main properties of the the LLC model, 
the MF approach, its temporal behaviour and bifurca¬ 
tion scenario. In Sec. m we define the system’s spatial 
geometry on a ring with classical nonlocal coupling and 
we numerically demonstrate the emergence of chimera 
states in dependence of the coupling range and the dis¬ 
tance from the critical point. In Sec. jlVj we introduce 
fractal connectivity matrices and give evidence of trav¬ 
eling chimeras. Additional results on connectivity with 
gaps are included in the Appendix. In the concluding 
section we summarize our results and discuss open prob¬ 
lems. 


II. PARAMETRIC STUDY OF THE LATTICE 
LIMIT CYCLE MODEL 

Originally, the LLC scheme [48] was devised as a highly 
nonlinear model for a cyclic reaction-diffusion process 
with predator-prey interactions among three particles A, 
Y, and S. The reaction scheme reads: 


2X + 2Y + S 

(la) 

X + S^2X 

(lb) 

Y + S^2S 

(Ic) 


where pi, p 2 , and ps are the interaction rates. 

The system Q being an ecological model was inspired 
by heterogeneous catalytic reactions, where single parti¬ 
cles are deposited on the sites of a catalyst surface. These 
particles interact with their neighbors or diffuse on the 
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surface. The catalyst is usually represented as a regular 
two-dimensional (2D) lattice (square, honeycomb, trian¬ 
gular), or even as a fractal lattice with impurities [58] . 
These catalytic applications indicate the need to repre¬ 
sent the empty lattice sites as a virtual species, which 
participates in the reaction and diffusion processes as 
well. The LLC model thus involves one virtual species, 
S, which represents the empty sites and two “normal” 
species, X and T, which are engaged in the three reac¬ 
tions. The scheme 0 describes that the interaction of 
two particles X and two particles Y turns one of the X- 
particles into Y while the other X-particle turns into S 
[Eq. ([la[)]. The positions (lattice sites) where the parti¬ 
cles are located may or may not be first neighbors, de¬ 
pending on the model (see Ref. [59]). Similarly, Eq. (lb) 
represents the birth of X if another X is found in an ad¬ 
jacent position, while Eq. (Ic) represents the cooperative 
death of a particle Y leaving an empty site S. 

This model has been extensively studied in recent 
years. Initially, it was implemented on a square lat¬ 
tice with single occupancy per lattice site using Kinetic 
Monte Carlo (KMC) simulations. Direct KMC realiza¬ 
tions on the 2D square lattice with nearest neighbor inter¬ 


actions produced intricate fractal patterns and local oscil¬ 
lations of the species concentrations [48] . Later on, long¬ 
distance diffusion was introduced as a mixing mechanism 
allowing the species to react with all particles within a 
specific range, thus giving them the possibility to change 
their places in the lattice at finite or infinite distances 
[59] . The model was also studied from the viewpoint of 
an abstract network of phases that was shown to have 
features of a scale-free network through calculations of 
the degree distribution and clustering coefficient m- 
In the ME approach, the system is described by three 
4th-order nonlinear differential equations for the tempo¬ 
ral total concentrations x, y, and s of the respective par¬ 
ticles X, K, and S. After applying the conservation con¬ 
dition X y -\- s = 1 , the LLC system is reduced to the 
following two equations: 


dx 

dt 

dt 


—2pix^y^ P 2 x{l — X — y) 
Pix'^y^ -P3y{l -x-y). 


( 2 a) 

( 2 b) 


The LLC system (Eq. 0 ) has four fixed points: one saddle-point Qi = ( 0 , 0 ), two other fixed points Q 2 = ( 0 , 1 ) 
with an unstable eigenvector direction along the ^-axis and Qs = (I 5 O) with a stable eigenvector direction along the 
x-axis (each having one zero eigenvalue [48]), and one nontrivial fixed point whose coordinates depend on the system 


parameters 




with K = y i+( 2 p 3 +P 2 )" 
V 27piP2P3 


In 


the parameter space (pi,P 2 ) for fixed ps, is either a stable node or a stable focus which becomes unstable through 
a supercritical Hopf bifurcation [48]. Because of the physical condition s > 0 the flow is always directed to the 
inside of the reaction simplex x> 0 ,p> 0 ,x + p<l which follows from mass-act ion kinetics. _ 


Hereafter, we fix the parameters p 2 = 0.5 and ps = 0.8, 
while Pi will be used as a control parameter determining 
the distance from the bifurcation point. Eor the above 
values of p 2 and ps, it can be found that the Hopf bifur¬ 
cation takes place at = 9.82. Eigure[^a) shows the 
limit cycle for three different values of pi > while 

the amplitude and period of the limit cycle as a function 
of Pi is plotted in Eigs. [^b) and (c), respectively. The 
X and y time series corresponding to the limit cycles of 
Eig. ^a) are depicted in Eigs. [^d) and (e), respectively. 
In the time series of the p-variable, in particular, the 
sharp spikes followed by long resting periods are visible 
for large values of pi. 

In the current study, we consider a setup based on ME 
oscillators arranged on a cyclic one-dimensional network 
(ring), interacting linearly with one another. The use of 
the ME dynamics on all network sites mimics the inter¬ 
actions between populations, where each population is 
composed of a mixture of x, p, and s. While the nodes 
interact with one another via nonlocal coupling as will 
be explained in the next section, the resident population 
follows the ME approach on each node given by Eqs. 0 . 
This is a different approach than the previous KMC simu¬ 


lations [36l SH] : where each node (lattice site) is occupied 
by a single particle. The properties of this particular ar¬ 
rangement of the populations is discussed in detail in the 
following section. 

III. MULTICHIMERA STATES IN THE 
LATTICE LIMIT CYCLE MODEL 

We consider N nonlocally coupled LLC oscillators on a 
one-dimensional ring according to the following scheme: 

dx (7 

= -2pixlyl + p2Xk{l - - 2/fc) + — ^ {xj - Xk) 

j=k — R 

(3a) 

, j=k+R 

=piXkyk-P3yk{i-Xk-yk) +^ 2 _^ [yj-yk)^ 

j=k — R 

(3b) 

where the subscript k refers to the node index, k = 
1,...,A, which has to be taken modulo N. Each os- 
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Figure 1: (Color online) (a) Limit cycle for various values 
of pi above the Hopf bifurcation. Fixed points Qi, Q 2 (not 
shown), Q 3 , and Q 4 (for pi = 10 ) are marked with crosses, 
(b) Peak-to-peak amplitude of the a:-variable and (c) Period 
of the limit cycle as a function of pi. (d) and (e) show the 
time series of x and y, respectively, corresponding to the limit 
cycles in (a). Other parameters: p 2 = 0.5 and ps = 0.8. 


cillator is coupled with its > 0 nearest neighbors on 
both sides with coupling strength a. This introduces 
nonlocality in the form of a ring topology as proposed in 
Ref. [ISIIGT]. There, the authors considered a rotational 
coupling matrix in order to achieve both direct and cross¬ 
coupling between the system variables. Our system is 
simpler in the sense that it contains only direct coupling 
in the species concentrations. This reflects mobility of 
the species between nodes and no cross-interactions of 
species between different nodes. 

A typical chimera state with two (in)coherent regions 
is shown in Fig.|^a), where the ^-concentration is plotted 
as snapshot. Initial conditions in all simulations herein 
are taken to be randomly distributed on an ellipsoid en¬ 
closed by the triangle with corners on the fixed points 
Qi, (52, and As a measure indicating the existence 
of a chimera state we employ the mean phase veloeity of 
each oscillator ujk = 27rM/e/AT, where is the number 
of periods of the kth oscillator during a time interval AT 
um- This quantity is depicted in Fig. |^b) and has 
the typical profile: flat, lower-valued in the coherent do¬ 
mains and arc-shaped, higher-valued, in the incoherent 
ones. Figure |^c) shows a snapshot of the (x, p)-plane 
(red dots) together with the limit cycle of the uncoupled 
system for the same parameter values. We observe that 
the coupling results in a limit cycle with smaller am¬ 
plitude than in the corresponding uncoupled system. Fi¬ 
nally, in the space-time plot of the x variable in Fig.j^d), 
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Figure 2 : (Color online) (a) Snapshot of the variable Xk of 
Eqs. § and (b) corresponding mean phase velocity profile, 
(c) Limit cycle of the uncoupled system (black solid line) and 
snapshot in the {x,y)-p\a,B.e (red dots), (d) Space-time plot 
of variable Xk- Parameters: pi = 300, p 2 = 0.5, ps = 0.8, 
N = 1000, R = 350, and cr = 0.015. 

it can be seen that the chimera state is fixed in space and 
has a period of approximately 18 time units, which corre¬ 
sponds to the period of the uncoupled system for a much 
lower value of pi « 40. 

The presence of chimera states in ecological and popu¬ 
lation systems can account for different dynamical states 
in different nodes of the interacting systems. For exam¬ 
ple, in population dynamics where communities exchange 
individuals, the density of inhabitants in some commu¬ 
nities may oscillate in phase, while in other communities 
it may be incoherent. Similarly, in computer networks, 
the density of users (or the data exchange) may oscillate 
coherently in some parts of the network, while in other 
parts it may behave asynchronously. 

In the following, we discuss how the coupling range 
R and the distance from the Hopf bifurcation, expressed 
through parameter pi , influence the form and multiplicity 
of chimera states. 


A. Impact of coupling range R 

Chimera states with multiple domains of incoherence 
and coherence have been reported in several works and 
are referred to as clustered chimera or multichimera 
states. It is known that they may be achieved through 
time delay [62] or by manipulating the range of the cou¬ 
pling between oscillators [Mni M- The range R of 
the coupling reflects the migration range of the different 
species in the system. Therefore, it is interesting to verify 
the formation of clustered chimera states in our system 
by varying this quantity. 

Figure shows typical chimera states and correspond- 
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Figure 3: (Color online) Snapshots of the x-variable (red 
dots) and corresponding ujk profiles (black dots) for various 
values of the coupling range: (a) R = 100, (b) R = 120, (c) 
R = 190, (d) R = 230, (e) R = 280, and (f) R = 410. Other 
parameters as in Fig. 


Figure 4: (Color online) Snapshots of the x-variable (red 
dots) and corresponding cuk profiles (black dots) at time t = 
3000, for a fixed coupling range R = 350 and various values 
of pn (a) pi = 20, (b) pi = 40, (c) pi = 90, (d) pi = 170, (e) 
pi = 210 and (f) pi = 370. Other parameters as in Fig. 


ing mean phase velocity profiles for increasing R. For 
low values of the coupling range, chimera states with 
up to 10 (in)coherent domains may be found [Fig. ia)]. 
For increasing R the multiplicity of the chimera state de¬ 
creases and, as R approaches the value of 280, a chimera 
state with 2 (in)coherent domains is found [Fig. [^e)]. 
This 2-chimera state persists over a large i7-interval up 
to ^ 400, where it starts to deform [Fig.j^f)]. The cor¬ 
responding mean phase velocity profiles follow the multi¬ 
plicity of the chimera states, exhibiting maximum values 
in the incoherent domains and constant values in the co¬ 
herent domains. 


B. Impact of bifurcation parameter pi 

As mentioned in the Introduction, for fixed rates p 2 
and p 3 , a Hopf bifurcation occurs at ^ 9.82 and a 
limit cycle starting from zero amplitude and finite pe¬ 
riod is born [see Figs. Bb) and (c)]. As pi increases, the 
limit cycle approaches the heteroclinic invariant mani¬ 
folds from Q 2 to Qi and from Qi to Qs which bound 
the basin of attraction of the limit cycle. The time to 
approach Qs becomes very large, and accordingly, its pe¬ 
riod diverges as pi increases further. Therefore, pi is a 
crucial parameter that introduces two time-scales in the 
system dynamics and determines the period of the indi¬ 
vidual oscillators. Physically, pi determines the rate of 
reactive process ( [Ta| ) and introduces a 4th-order term in 
the equations, since it requires simultaneous interaction 
of four species (particles). This parameter takes high val¬ 
ues compared to the other two, p 2 and ps, to compensate 
the contribution of the term x^y‘^. 

Figure shows typical snapshots of the x-variable and 


corresponding mean phase velocity profiles for fixed cou¬ 
pling range R = 350 and increasing pi. For values of pi 
close to the bifurcation point [Fig.j^a)], no chimera state 
is observed but, instead, a mixed state with no obvious 
spatial structure is present where some oscillators clus¬ 
ter, while others move incoherently. This is also reflected 
in the mean phase velocity profile depicted in the lower 
panel of Fig. [^a), which does not depict the character¬ 
istic arc of a chimera state. The amplitude of the oscil¬ 
lations for this parameter values is small. As we move 
away from the bifurcation point, the oscillators are or¬ 
ganized in one coherent and one incoherent region, thus, 
forming a 1-chimera state [Figs. Qb) and (c)]. The ujk- 
profile starts to obtain its typical shape, which becomes 
more prominent as pi attains higher values. Then, a 2- 
chimera state appears around pi = 190 [Figs. ®e) and 

(f)]- 

A careful look at Figs, [^c)-(f) reveals the scenario of 
how a chimera state with two incoherent parts is formed. 
The mean phase velocity profile develops a small dip, 
which increases and eventually reaches the level of the co¬ 
herent part as the parameter pi increases. Similar scenar¬ 
ios have been observed for nonlocally coupled FitzHugh- 
Nagumo systems m- 

The impact of pi is also depicted in Fig. where the 
space-time plots for the variable x (top panel) and the 
corresponding local order parameter (lower panel) are 
shown for fixed R = 350. The values of pi are chosen as 
Pi = 90,170,210, and 370 in panels [^a) to[^d), respec¬ 
tively, which corresponds to the snapshots in Figs , [^c) - 
(f). The local order parameter is defined as follows IbTTbdj : 


= 


1 


\j-k\<5 


k = l,...,N, 


(4) 
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Figure 5: (Color online) Space-time plots of the x-variable (top) and order parameter (bottom) for fixed R = 350 and increasing 
Pi- (a) Pi = 90, (b) Pi = 170, (c) pi = 210, and (d) pi = 370. Other parameters as in Fig.[^ 


where Qj = arctan[(^j — yQ^)/{xj — xq^)] denotes the 
geometric phase of the jth LLC unit and {xq^^yq^) are 
the coordinates of the nontrivial fixed point Q/i of the 
uncoupled system. We use a spatial average with a win¬ 
dow size of ^ = 25 elements. The local order parameter 
Zk close to unity indicates that the kth unit belongs to 
the coherent part of the chimera state, while is less 
than 1 for incoherent parts. 

The lower panels of Fig. depict the local order pa¬ 
rameter in the time interval t G [2700, 3000], where bright 
(yellow) color denotes the coherent regions. A station¬ 
ary 2-chimera state is obtained for high values of pi, as 
shown in Fig. [^d). As one moves closer to the Hopf bi¬ 
furcation point, the two incoherent regions merge into a 
large one, while the size (number of oscillators) of the in¬ 
coherent region increases at the expense of the coherent 
ones [Fig. I^a) and (b)]. 


with two gaps. The introduction of regular gaps induces 
merging and splitting of the (in) coherent parts of the 
chimera state but does not produce other qualitatively 
different features. Representative cases are shown in the 
Appendix 


(a) 


n=0 


n=l 


n=2 




IV. TRAVELING CHIMERAS FOR 
HIERARCHICAL CONNECTIVITY 

In a previous study on the FitzHugh-Nagumo model 
we have demonstrated that if the connectivity matrices 
have hierarchical form, the chimeras also occur as nested 
structures m- Here, in addition to these structures, we 
show evidence of a different phenomenon related to this 
connectivity, namely, the occurrence of traveling coher¬ 
ent /incoherent regions. 

The hierarchical coupling structure involves connec¬ 
tivity gaps. The simplest example is a coupling matrix 


Figure 6: (Color online) Schematic representation of the 

hierarchical connectivity using the triadic Cantor set. (a) 
“bottom-up” construction of the triadic Cantor set, (b) linear 
connectivity arrangement and (c) ring connectivity arrange¬ 
ment for one reference node. 

Before we discuss the effect of the hierarchical connec¬ 
tivity, we summarize the procedure to generate the corre¬ 
sponding network. Figure illustrates this procedure in 
a schematic diagram [65] . Let us consider a connectivity 
matrix produced by an initiation string or base of size 5, 
e.g. 5 = 3 with a base pattern S =“101” in Fig. [^ con- 
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taining ci times the symbol 1 and Cq = b — Ci times the 
symbol 0. The specific arrangement of the symbols on 
the initiation string is essential, because the hierarchical 
connectivity pattern is constructed as the nth iteration 
of the initiation string replacing 1 by the base pattern S 
and 0 by 6 times the symbol 0. Thus, the length of the 
string defining the connectivity of the network is equal to 
the system size N = h'^ and contains times the symbol 
1 and N — the symbol 0. The limiting set, which is 
produced when the number of iterations n ^ oc, is a frac¬ 
tal Cantor set and has fractal dimension df = lnci/ln6. 
Since d/ is formally defined for infinite systems only, the 
finite size sets used here are called hierarchical, because 
they have been constructed based on a hierarchical algo¬ 
rithm. As a last step, we obtain the total connectivity 
of the network by index shift. As we shift the iterated 
string along the ring, we can determine the connectivity 
for each node. In this construction, each node has pre¬ 
cisely Cl links to other nodes hierarchically arranged in 
a unique pattern and directed. The corresponding con¬ 
nectivity matrix is given by 

{ 1 if both nodes k and I belong to the 
Cantor set obtained by n iterations 
0 elsewhere. 

(5) 

This connectivity matrix is of size x b'^ and contains 
a hierarchical distribution of gaps with a variety of sizes. 

In the following, we use the base size 6 = 6 and n = 4 
iteration steps producing a system of size N = 6^ = 1296. 
The considered initiation strings are S =“001111” and 
S =“110111”, which iterated n = 4 times produce the 
connectivity patterns, consisting of 256 and 625 times 
the symbol 1 and 1040 and 671 times the symbol 0, re¬ 
spectively. 
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Figure 7: (Color online) Space-time plots of two traveling 
chimeras for (a) df = In 4/In 6 = 0.774 and initiation string 
“001111” and (b) df = In5/ln6 = 0.898 and initiation string 
“110111”. All other parameters as in Fig. 


position are presented. The space-time plots in Fig. 
also reveal an additional internal nested structure in both 
cases which are the imprints of the hierarchical connectiv¬ 
ity. Similar nested chimera states have been previously 
demonstrated by the FitzHugh-Nagumo system without 
the additional feature of traveling (in)coherent parts [64] . 




Figure 8: (Color online) (a) Position of the oscillator display¬ 
ing maximum x-value with time for a = 0.015 and (b) Fourier 
transforms for cr = 0.015 (red spectrum) and a = 0.020 (black 
spectrum). The inset in b) shows in detail the lowest part 
of the two spectra where the frequencies associated with the 
traveling motion around the ring are discerned. All other 
parameters as in Fig.[^b). 

In the case of traveling chimeras, the mean phase ve¬ 
locity is not a good measure for coherence, because each 
oscillator spends part of the time in the coherent regions 
and part of the time in the incoherent ones. An alter¬ 
native way to distinguish between stationary and travel¬ 
ing chimeras is to locate the nodes attaining maximum 
and minimum x-values (or ^-values) at each time step. 
Given that the incoherent regions are not localized but 
travel along the ring in time, we expect that the node 
(position on the ring) which achieves maximum x-values 
(or ^-values) is also subject to change. In Fig. |^a), we 
record the node, JmaxCO, that exhibits maximum x-value 
at time t, which yields the condition: 


The motion of the (in) coherent regions in traveling 
chimeras are depicted in the space-time plots. Fig. Ilia) 
and (b). From this figure we can observe a continuous 
translation of the (in)coherent parts along the ring, a 
motion which is not observed, for example, in Fig. 
where exemplary space-time plots of chimeras with fixed 


Similarly, one can define the element which ex¬ 

hibits minimum value at time t. The position Jmax(^) (or 
^min(^)) shows a pcriodicity, since each element succes¬ 
sively passes through the state of maximum amplitude. 





































































The period T needed for a single element to attain the 
state of maximum amplitude can be calculated using the 
Fourier transform of the time series Jmax (^)5 Fig- |^^)- 
The Fourier transform, Fig. 8(b), presents two distinc¬ 
tive maxima. For example, for a = 0.015 we first observe 
in Fig. 8(b) (red line) one maximum of high frequency 
/osc = 0.056, and short period T^sc = l//osc = 17.86 
which correspond to the average frequency of the oscilla¬ 
tors. By looking at the low frequencies (inset) we observe 
an additional maximum at frequency /tr = 0.00071 with 
long period Ttr = 1408, which is associated with the trav¬ 
eling motion around the ring. Note that Tqsc is smaller 
than the period of the uncoupled oscillator T(pi = 300) 
[cf. Fig. [^c)], because each unit is - from time to time 
- part of the incoherent domain, which oscillates faster 
than the coherent region. Furthermore, the traveling pe¬ 
riod Ttr can also be inferred from the time series Jmax 
shown in Fig. |^a). The black curves in Figure [^b) re¬ 
fer to a larger coupling strength a = 0.02 than the one 
used in all previous plots (a = 0.015). Comparing the 
spectra for the two a values in Fig.j^b) we find that the 
frequency ftr moves to lower values as a increases. 

This is further detailed in Fig. [^that shows the trav¬ 
eling speed Vtj. of the domains. During one period Ttr, 
the coherent and incoherent regions cover the entire ring 
consisting of N elements, thus, 

t'tr = N/Ttr = Nftr (7) 

The chimera traveling speed Vtr decreases for increas¬ 
ing coupling strength cr, i.e. as the coupling con¬ 
stant becomes larger, the incoherent regions slow down 
around the ring. We recall here the results of Sec. |III| 
where chimera states are only observed for low coupling 
strengths in the LLC model. 



Figure 9: (Color online) Traveling speed utr of (in)coherent 
regions around the ring. All other parameters as in Fig.j^b). 

As alternative measure for the calculation of the trav¬ 
eling speed we may use the spectrum of single oscillators 
(not shown). Since each oscillator belongs periodically 
either to the coherent or the incoherent regions, it will 
participate once during each period in the regions which 
demonstrate maximum (or minimum) amplitude. Then 


by calculating the spectrum of a single oscillator and se¬ 
lecting the lowest frequency we obtain an estimation of 
the chimera traveling speed based on one oscillator. 
Average over the N elements needs to be taken to repre¬ 
sent the overall chimera speed Vtr- 


V. CONCLUSIONS 

In the present work, we have demonstrated for the first 
time the existence of chimera states in a population dy¬ 
namics model, reflecting the various tendencies of local 
communities to behave (oscillate) coherently or incoher¬ 
ently within the same network. We have shown that the 
dynamics of the nonlocally coupled system strongly de¬ 
pends on the character of the local behavior of the nodes. 
When the parameters of the individual LLC nodes are 
close to the Hopf bifurcation, chimera states cannot be 
observed in the system. In this regime, all strongly non¬ 
linear 4th-order terms become weak and the correspond¬ 
ing limit cycle oscillations lose their spiking form resem¬ 
bling simple harmonic motion of small amplitude. When 
the internal bifurcation parameter of each node increases 
away from the Hopf bifurcation, the trajectories approach 
the heteroclinic invariant manifolds of the saddle points 
producing spikes followed by long resting periods. This 
type of local dynamics allows for the existence of a va¬ 
riety of chimera and multichimera states, which depend 
both on the parameters of the individual system, and the 
coupling. 

While the introduction of single gaps in the con¬ 
nectivity matrix induces merging and splitting of the 
(in)coherent parts of the chimera states (see Ap¬ 
pendix 0 , hierarchical arrangement of gaps results in 
the emergence of nontrivial phenomena in which the 
(in)coherent regions show nested structures and travel 
along the ring, keeping their profiles statistically stable in 
time. Moreover, we have found that the speed of this mo¬ 
tion decreases with increasing coupling strength. Com¬ 
plex nested chimera structures, when regarded from the 
viewpoint of population dynamics, show the rich organi¬ 
zation which can emerge in communities of nonlinearly 
interacting populations. 
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Appendix A: Chimera states for connectivity 
matrices with gaps 


Previous studies have demonstrated that the intro¬ 
duction of gaps in a system of coupled elements mod¬ 
ifies the chimera states, producing different coherence- 
incoherence patterns [64] . We test this hypothesis in the 
case of the LLC model and investigate whether this ef¬ 
fect is generic or exclusive to specific models. For this 
we consider an adjacency matrix, {Cki}, k^l = 1 ,...,A, 
with two gaps, one to the left and one to the right of 
the node k, around which the connectivity is described, 
called reference node. Generally, the size of the gaps and 
the sizes of the linked regions can vary independently on 
the left and on the right of the reference node. The form 
of matrix element Cki is: 


Cki 


1 if k — I <C k Ri 

OT k Ri + Gr I k Ri + Gr + R 2 
OT k — Rs — Gr — R4<l<k — Rs — Gr 
0 elsewhere, 

(Al) 


where all indices are taken modulo N. This general form 
indicates that the node k is linked to two groups of nodes 
to the right with sizes Ri and R2 which are separated 
by a gap of size Gr. A similar connectivity arrange¬ 
ment with linked groups of sizes Rs and R4 and gap size 
Gr applies to the left of each node. Symmetr ic co nnec- 
tivity is realized as a particular case of Eq. with 

i ?3 = i?i, Gr = Gr and i ?4 = i? 2 - The connectivity 
matrix Eq. can be generalized further, to present 

multiple gaps to the left and to the right interrupting the 
linked elements. An asymmetric connectivity matrix, for 
instance, describes a directed network. 

We now test the influence of the gap size using asym¬ 
metric connectivity matrices with links extending only 
to the right of each node. Because symmetry in the con¬ 
nectivity does not influence the relative fraction of nodes 
belonging to the (in)coherent part of the chimera, we 
only consider links to the right part of each node. This 
reduction is introduced to keep the number of parameters 
minimal; by using i ?3 = i ?4 = Gl = 0 we only vary three 
parameters related to connectivity, i?i, Gr and R 2 . 

To investigate the influence of the gap size on coher¬ 
ence, we keep the sizes of the two linked regions Ri and 
R 2 constant, and vary only the gap size Gr in the inter¬ 
val 0 < Gr < i?i + i? 2 - In Eig. we plot the typical 
chimera profiles as we increase the gap size Gr. The sys¬ 
tem size is kept to N = 1000 oscillators, with parameter 
values Pi = 300, p 2 = 0.5, p 3 = 0.8, a = 0.015. The total 
number of connections is also kept fixed, R = 400, and is 
divided into two regions Ri = 100 and R 2 = 300. 

In Eig. we observe changes in the chimera multi¬ 
plicity, ranging from one to three (in)coherent regions, 
merging and splitting of coherent and incoherent regions 
and shifting of their position in space. The change of the 
chimera profiles can be quantified statistically by calcu¬ 
lating the following measures of coherence m- a) the av- 




Figure 10: (Color online) Chimera states (left) and cor¬ 

responding mean phase velocity profiles (right) for gap size 
variation. The gap sizes are indicated in the right panels, 
while Ri = 100 and R 2 = 300. All simulations start from the 
same initial conditions. Other parameters as in Fig. 


erage mean phase velocity of the coherent parts, (cjcoh), 
b) the maximum difference of the mean phase velocities 
A(j = uJmax — ^mini c) the fraction of oscillators belong¬ 
ing to the incoherent parts, A^incoh, and d) the extensive 
measure of incoherence Mincoh- The relative size Wncoh 
of the incoherent parts of the chimera state is calculated 
as: 


-^incoh — ^ ^ (^coh) c) 

^ k=l 

where 0 is the step function which takes the value 1 when 
its argument takes positive values and zero otherwise, c 
is a small tolerance, which in this case we set to 0.05. We 
also define the extensive, cumulative size Mincoh of tho 
incoherent parts as: 


N 

-^incoh — ^ ^ \ {^k (^coh))| 

k=l 

This is an extensive measure which represents the area 
below the arcs in the mean phase velocity profiles. Mincoh 
is a measure of incoherence and is equal to zero for fully 
coherent states. 

In Fig. the four coherence measures a plotted, av¬ 
eraged over 10 different initial conditions. Variation of 
the gap size shows a slight increase in the ujcoh value for 
small values of G, which soon attains a constant value, 
independent of the gap size. The other three measures 
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Figure 11: (Color online) The four measures of coherence as 
a function of the gap size G: a) (o^coh), b) Alu = comax 
c) Mncoh, d) Mincoh- Avcragcs are taken over 10 different sets 
of initial conditions. Parameter values as in Fig. |10[ 


of coherence do not show any appreciable change with 
variations of the gap size. Especially the number of inco¬ 
herent and coherent oscillators A^incoh and N — A^incoh do 
not change drastically with G, although the number of 
coherent and incoherent regions may change, as indicated 
by Fig. This means that when the gap size changes 
there is a redistribution of the incoherent oscillators in 


one large region (Fig. [1Q|^c)) or in two or three smaller 
ones (Fig. [^a),(b),(d)y7so that the total number of in¬ 
coherent oscillators is conserved (see Fig. 11). 


Earlier studies for connectivity matrices with gaps in 
the EitzHugh-Nagumo system have demonstrated that 
the position of the gap is important and that it affects the 
chimera properties modifying all four measures of coher¬ 
ence M- To test this for the LLC model we produce the 
snapshots and corresponding mean phase velocity profiles 
varying the sizes of the connectivity regions Ri and R 2 
but keeping the total number of links = 400 fixed, 

while introducing a gap of constant size Gr = G = 100 
at various positions between the links. The results are 
depicted in Eigurep!^ 


In contrast to the EitzHugh Nagumo system, here we 
do not observe a systematic dependence in the number 
of the (in)coherent regions as the gap of constant size 
moves away from the reference node. To further verify 
this, we plot in Eig. the four measures of coherence 
for different values of Ri. Averages are taken over 10 
different different initial conditions. All four measures 
seem to be constant (up to fluctuations) and they do not 
show any systematic change as a function of Ri. 
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Figure 12 : (Color online) Chimera states and corresponding 
mean phase velocity with variations on the position of a gap 
of constant size Gr — 100. (a) Ri — 50,i?2 = 350, (b) 

Ri — 200, = 200, (c) = 300, i ?2 = 100, and d) Ri — 

400, i ?2 = 0. All runs start from the same initial conditions. 
Other parameters as in Fig. 



Figure 13: (Color online) The four measures of coherence as a 
function of the inner connectivity radius Ri , while Ri-\-R 2 = 
400: (a) (cJcoh)? (b) — ^max ^min-j (c) Aincoh 5 (d) Afincoh • 

The gap size is kept to Gr = 100. Other parameter values 
as in Fig.[^ Averages are calculated over 10 different initial 
conditions. 





























